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dynamics  and  intramolecular  energy  transfer  processes  in  larger  polyatomic 
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elementary  reactions  of  importance  in  the  chemical  vapor  deposition  of  silicon  from 
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decomposition  reactions  of  Siha,  SiH4,  Si2H4,  and  SiaHe,  a  variety  of  nucleation 
reactions  leading  to  silicon  clusters  such  as  Sia,  Sis,  and  Si4,  and  polymerization 
reactions  leading  to  the  higher  silanes  and  silenes.  Heterogeneous  processes 
studied  include  sticking  probabilities,  scattering,  and  chemical  reactions  of  Ha  and 
SiHa  on  Si(IOO),  Si(111)'(1x1),  and  reconstructed  Si(111)-{7x7)  surfaces.  Surface 
diffusion  of  hydrogen  and  silicon  atoms  on  Si{111)-(1x1)  and  Si(111)-(7x7)  have  also 
received  considerable  attention. 

In  addition  to  the  investigation  of  the  dynamics  of  the  above  systems,  we  have  also 
developed  several  new  theoretical  methods  required  to  study  such  complex  systems. 
These  methods  include;  (1)  methods  for  obtaining  sufficiently  accurate  global 
potential-energy  surfaces  for  use  in  trajectory  calculations  and  in  variational  phase- 
space  studies,  (2)  velocity  reset  procedures  to  simulate  the  effects  of  relaxation  to  bulk 
phonon  modes  of  a  surface,  (3)  perturbation  methods  that  permit  accurate  gas-surface 
scattering  studies  to  be  carried  out  in  a  significantly  reduced  dimensional  framework, 
(4)  semi-empirical  procedures  for  the  computation  of  two-dimensional  surface 
tunneling  rates  that  take  into  account  the  effect  of  all  of  the  surface  phonon  modes  and 
possib'e  tunneling  paths,  (5)  h'ghly  efficient  variational  phase-space  methods  that 
incorporate  exact  evaluation  of  the  integrals  over  momentum  space  and  Efficient 
Microcanonical  Sampling  techniques,  and  (6)  a  new  projection  method  for  the 
computation  of  intramolecular  energy  transfer  rates  that  obviates  the  necessity  to 
arbitrarily  define  a  "bond"  or  "mode  energy".  These  methods  have  been  tested  by  the 
silicon  research,  by  studies  of  scattering  of  NO  from  Ag(111)  surfaces,  and  by  the 
investigation  of  the  bimolecular  reaction  mechanism  for  CaH4  +  Fa  and  the 
unimolecular  reaction  of  1 ,2-difluoroethane.  As  a  result  of  these  investigations,  we 
have  discoved  the  presence  of  striking  non-statistical  effects  in  several  gas-phase, 
unimolecular  reactions  that  call  into  question  the  general  use  of  statistical  theories  to 
treat  such  reactions.  For  example,  we  have  found  that  simple  bond-fission  reactions  of 
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Abstract: 


This  Final  Report  for  Grant  AFOSR-89-0085  summarizes  the  results  of 
research  conducted  under  AFOSR  support  with  particular  emphasis  on  the 
investigations  carried  out  during  the  three-year  grant  period  for  AFOSR-89- 
0085.  The  research  reviewed  includes  homogeneous  and  heterogeneous 
processes  of  particular  importance  in  the  chemical  vapor  deposition  (CVD)  of 
silicon  from  silane  and  disilane,  the  study  of  chemical  processes  occurring 
under  conditions  of  close  confinement,  non-statistical  dynamics  and 
intramolecular  energy  transfer  processes.  New  methods  for  (1)  obtaining 
potential  energy  surfaces  for  highly  complex  systems,  (2)  simulation  of  the 
effects  of  relaxation  to  the  bulk  in  surface  systems,  (3)  perturbation  studies  of 
gas-surface  scattering,  (4)  computation  of  two-dimensional  surface  tunneling 
rates,  (5)  highly  efficient  variational  phase  space  theory  calculation  of 
microconical  unimolecular  reaction  rates,  and  (6)  the  computation  of 
intramolecular  vibrational  relaxation  rates  are  also  described. 
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I,  Introduction: 


This  report  summarizes  the  research  programs  conducted  with  AFOSR 
support  under  Grant  AFOSR-89-0085.  For  purposes  of  clarity,  the  research  is 
discussed  in  the  context  of  work  completed  under  previous  AFOSR  support. 

The  program  involves  primarily  the  theoretical  study  of  the  elementary 
reactions  of  importance  in  the  chemical  vapor  deposition  of  silicon  from  silane. 
This  includes  the  investigation  of  homogeneous  unimolecular  decomposition 
riactions  of  SiHa,  SiH4,  Si2H4,  and  Si2H6,  a  variety  of  nucleation  reactions 
leading  to  silicon  clusters  such  as  Si2.  Sia,  and  Si4,  and  polymerization 
reactions  leading  to  the  higher  silanes  and  silenes.  Heterogeneous  processes 
under  study  include  sticking  probabilities,  scattering,  and  chemical  reactions  of 
H2  and  SiH2  on  Si(100),  Si(1 1 1)-(1x1),  and  reconstructed  Si(1 1 1)-(7x7) 
surfaces.  Surface  diffusion  of  hydrogen  and  silicon  atoms  on  Si(111)-(1x1)  and 
Si(1 1 1)-(7x7)  have  also  received  considerable  attention. 

In  addition  to  the  investigation  of  the  dynamics  of  the  above  systems,  we 
have  also  developed  several  new  theoretical  methods  required  to  study  such 
complex  systems.  These  methods  include:  (1)  methods  for  obtaining  sufficiently 
accurate  global  potential-energy  surfaces  for  use  in  trajectory  calculations  and 
in  variational  phase-space  studies,  (2)  velocity  reset  procedures  to  simulate  the 
effects  of  relaxation  to  bulk  phonon  modes  of  a  surface,  (3)  perturbation 
methods  that  permit  accurate  gas-surface  scattering  studies  to  be  carried  out  in 
a  significantly  reduced  dimensional  framework,  (4)  semi-empirical  procedures 
for  the  computation  of  two-dimensional  surface  tunneling  rates  that  take  into 
account  the  effect  of  all  of  the  surface  phonon  modes  and  possible  tunneling 
paths,  (5)  highly  efficient  variational  phase-space  methods  that  incorporate 
exact  evaluation  of  the  integrals  over  momentum  space  and  Efficient 
Microcanonical  Sampling  techniques,  and  (6)  a  new  projection  method  for  the 
computation  of  intramolecular  energy  transfer  rates  that  obviates  the  necessity 
to  arbitrarily  define  a  "bond"  or  "mode  energy".  These  methods  have  been 
tested  by  the  silicon  research,  by  studies  of  scattering  of  NO  from  Ag(111) 
surfaces,  and  by  the  investigation  of  the  bimolecular  reaction  mechanism  for 
C2H4  +  F2  and  the  unimolecular  reaction  of  1 ,2-difluoroethane.  As  a  result  of 
these  investigations,  we  have  discoved  the  presence  of  striking  non-statistical 
effects  in  several  gas-phase,  unimolecular  reactions  that  call  into  question  the 
general  use  of  statistical  theories  to  treat  such  reactions. 


In  Section  II,  we  review  the  previous  and  current  progress  of  our 
research  program.  The  section  concludes  with  a  listing  of  all  papers  and 
presentations  made  under  AFOSR  sponsorship.  The  personnel  associated  with 
the  program  are  described  in  Section  III. 

For  the  reader's  convenience,  we  have  concluded  most  sections  and 
subsections  with  a  brief  summary.  These  summaries  may  be  used  to  obtain  a 
rapid  overview  of  the  entire  proposal. 


II.  Review  of  Previous  Work 


IIA.  Overview 

We  have  been  conducting  theoretical  studies  of  the  elementary 
homogeneous  and  heterogeneous  processes  involved  in  the  chemical  vapor 
deposition  (CVD)  of  silicon  from  silane  and  disilane.  The  homogeneous  gas- 
phase  reactions  investigated  to  date  include  the  key  unimolecular  dissociation 
processes  of  SiaHe,  Si2H4,  SiH4,  and  SiHa.  In  each  case,  the  calculations 
have  treated  all  energetically  accessible  decomposition  channels.  Similar 
studies  have  also  been  carried  out  on  1 ,2'difluoroethane  as  part  of  a  program 
designed  to  study  reactions  occurring  under  highly  confined  conditions. 
Bimolecular  reactions  that  constitute  important  steps  in  polymerization 
processes  leading  to  higher  silanes  have  received  considerable  attention. 
Silicon  clustering  reactions  that  lead  to  the  formation  of  Si2,  Sia,  and  Si4 
aggregates  have  been  examined  and  their  relative  importance  in  promoting 
silicon  crystal  growth  in  the  CVD  experiments  determined. 

The  heterogeneous  silicon  surface  processes  that  have  been 
investigated  to  date  include  studies  of  sticking,  dissociative  chemisorption, 
energy  transfer,  surface  diffusion,  surface  relaxation  to  bulk  modes,  chemical 
reactions,  and  tunneling  effects  of  H,  H2,  Si,  Si2,  Sia,  Si4,  and  SiHa  species  on 
Si(100),  Si(111)-(1x1),  and  reconstructed  Si(111H7x7)  surfaces. 

As  mentioned  above,  we  have  initiated  investigations  of  chemical 
processes  that  occur  in  very  confined  spaces.  The  study  of  such  reactions  has 
been  listed  as  one  of  the  key  problems  in  silicon  chemistry.^  In  order  to 
determine  the  effects  of  close  confinement  upon  the  reactions  in  question,  we 
have  first  conducted  detailed  dynamical  studies  of  the  corresponding  gas- 
phase  reactions  to  provide  a  basis  for  comparison  with  the  results  obtained 
when  the  reactions  occur  in  a  confined  environment.  The  gas-phase  dynamics 
predicted  by  these  investigations  have  also  provided  us  with  a  means  of 
evaluating  the  accuracy  of  our  methods  for  the  development  of  global  potential- 
energy  surfaces. 

As  a  result  of  our  extensive  investigations  of  microcanonical 
unimolecular  decomposition  reactions,  we  have  recently  discovered  the 
existence  of  non-statistical  effects  that  call  into  question  the  use  of  transition- 
state  methods  for  the  study  of  such  reactions.  To  date,  our  investigations, 


coupled  with  some  recently  reported  experimental  results,  suggest  that  these 
non-statistical  effects  may  be  widespread  in  reactions  of  polyatomic  molecules. 
By  comparing  the  detailed  molecular  and  intramolecular  energy  transfer 
dynamics  of  the  2-chloroethyl  radical  (which  behaves  statistically)  with  those  for 
1 ,2-d'f!unroethane  and  disilane  (which  are  nonstatistical),  we  have  been  able  to 
formuiaia  ome  general  principles  that  permit  perdictions  of  nonstatistical 
dynamics  to  be  made  in  advance  of  the  experiments  and/or  calculations. 

The  chemistry  involved  in  the  above  systems  is,  in  many  cases,  very 
complex.  Often  several  competing  reaction  channels  must  be  considered.  This 
has  required  the  use  of  a  wide  variety  of  theoretical  methods  that  include 
classical  trajectories^,  variational  phase-space  theory  methods  such  as  those 
developed  by  Truhlar,  Garrett,  and  coworkers^  and  by  Voter  ana  DolK,  Monte 
Carlo  methods  with  efficient  microcanonical  sampling  procedures  (EMS)®-^, 
wavepacket  calculations,  perturbation  methods®'^,  projection  procedures'* o  to 
follow  intramolecular  energy  transfer  processes,  Langevin  methods  developed 
by  Aldeman,  Doll,  and  coworkers^*  and  by  Lucchese  and  Tully*^,  and  velocity 
reset  methods  such  as  those  reported  by  Riley,  Coltrin,  and  Diestler^®.  In 
addition,  the  need  to  obtain  sufficiently  accurate  potential-energy  surfaces  has 
required  the  development  of  procedures  that  permit  the  results  of  ab  initio  SCF, 
Cl.  and  Mbiler-Plesset  fourth-order  perturbation  calculations  (MP4)  to  be 
combined  with  pertinent  experimental  data  to  produce  the  required  global 
potential-energy  surfaces.  In  order  to  study  surface  tunneling  effects,  we  have 
developed  procedures  that  permit  reaction-path  methods  to  be  coupled  with 
variational  phase-space  theory*'^. 

The  general  thrust  of  our  research  has  been  to  determine  the  extent  to 
which  existing  methods  are  useful  in  the  examination  of  complex  chemical 
systems  and,  where  necessary,  to  develop  suitable  extensions  of  existing 
methods  or  entirely  new  procedures.  In  this  section,  we  review  our 
accomplishments  to  date.  The  new  methods  developed  and  the  principal 
results  of  each  investigation  are  briefly  summarized.  At  the  end  of  the  section, 
we  give  a  complete  listing  of  all  papers  published  and  oral  presentations  given. 

ILfii _ Potential-Energy  Hvpersurfaces 

The  most  crucial  and  difficult  part  of  any  dynamics  study  for  a  complex 
system  is  the  development  of  a  potential-energy  hypersurface  that  is  sufficiently 


close  to  the  potential  for  the  experimental  system  under  investigation  that  the 
results  of  the  dynamics  calculations  may  be  regarded  as  meaningful.  Such 
surfaces  must  accurately  represent  all  open  reaction  channels.  Since  it  is 
obviously  impossible  to  carry  out  converged,  multi-referenced,  complete 
configuration  interaction  calculations  at  a  sufficiently  large  number  of 
configuration  points  to  completely  characterize  the  hypersurface,  the  problem 
can  never  be  "solved".  However,  our  experience  has  shown  that  this  level  of 
accuracy  is  not  necessary  to  obtain  useful  results.  If  cetrain  key  elements  of  the 
potential  surface  are  accurately  described  and  if  the  choice  of  functional  forms 
used  to  the  fit  the  available  experimental  and  ab  initio  data  is  based  on  well- 
judged  chemical  and  physical  considerations,  then  the  computed  dynamics, 
rates,  and  mechanisms  for  the  system  are  often  in  good  to  excellent  accord  with 
experiment  even  though  many  of  the  finer  details  of  the  surface  may  be 
incorrect. 

The  minimum  data  base  required  for  the  formulation  of  a  global  potential 
hypersurface  contains  data  on  the  thermochemistry  for  all  important  reactant 
channels,  fundamental  vibrational  frequencies  and  equilibrium  geometric 
structures  for  ail  reactants  and  products,  and  potential-energy  barriers  for  the 
reactions  of  importance.  It  is  also  very  useful  to  have  data  concerning  the 
structure  and  vibrational  frequencies  of  the  transition  states.  Much  of  this  data  is 
usually  available  from  experimental  measurement.  Data  related  to  the 
transition  states  and  radical  species  can  be  obtained  from  electronic  structure 
calculations. 

Once  the  required  data  base  is  in-hand,  a  suitable  analytical  function  for 
the  system  potential  that  can  be  accurately  fitted  to  the  data  must  be  developed. 
In  principle,  the  total  potential  may  be  written  as  a  many-body  expansion  of  the 
form 

(1)  V  =  Z  Vi, 

where  the  Vj  are  i-body  terms  and  V  includes  the  sum  of  all  possible  such  terms. 
Murrell  and  co-workers^5  have  made  effective  use  of  such  an  expansion  for 
systems  in  which  the  number  of  atoms,  N,  is  <  4.  For  more  complex  systems, 
however,  Eq.(1 )  becomes  computationally  intractable.  In  fact,  even  for  the  case 
N  =  4,  it  is  extremely  difficult  to  use  Eq.(1).  For  systems  as  complex  as  those  of 
interest  here,  it  is  essentially  impossible. 


During  the  past  three  years,  we  have  developed  an  approach  to  these 
problems  that  has  proven  itself  to  be  very  satisfactory.  Basically,  we  employ 
many-body,  parametrized  functional  forms  suggested  by  chemical  and  physical 
considerations.  This  analytical  potential  is  then  fitted  to  the  data  base  using  an 
iterative  scheme  involving  simultaneous  parameter  adjustment.  Experience 
has  shown  that  excellent  fits  can  be  achieved  using  such  methods. 
Furthermore,  we  also  find  that  the  method  often  produces  acceptable  accuracy 
with  regard  to  topographical  features  of  the  surface  not  directly  addressed  in  the 
fitting  process.  Numerous  direct  comparisons  of  the  results  of  dynamics 
calculations  on  potential  surfaces  obtained  in  this  manner  with  experimental 
data  have  shown  the  extent  of  agreement  to  be  generally  good,  and  often 
excellent.  Several  examples  of  such  agreement  are  given  in  the  following 
sections.  Although  tedious  and  often  difficult  to  implement,  we  believe  that  such 
a  procedure  is  the  only  currently  available  route  to  ensure  that  the  potential  for 
the  model  system  will  be  sufficiently  close  to  that  for  the  actual  experimental 
system  to  permit  meaningful  dynamics  calculations  to  be  carried  out. 

We  have  previously  applied  the  above  method  in  the  development  of 
potential  hypersurfaces  for  SiHa^®,  SU^®,  Si2H4^®,  SiaHe^®,  CaH420, 
and  HaFC-CHaF20.2i.  in  each  case,  these  surfaces  correctly  predict  a  wide 
array  of  structural,  energetic,  and  dynamical  data  for  the  system.  The  results 
discussed  in  the  following  sections  serve  as  illustrations  of  the  level  of  accuracy 
that  we  have  been  able  to  obtain  using  these  methods. 

S-umaary: 

Global  potential-energy  hypersurfaces  are  obtained  by 
employing  many-body,  parametrized  functional  forms  suggested  by 
chemical  and  physical  considerations  to  fit  a  data  base  consisting 
of  energies,  structures,  and  potentiai  surface  curvatures  for 
reactants,  products,  and  transition  states  obtained  from 
experimentai  measurements  and  the  resuits  of  ab  initio  eiectronic 
structure  caiculations.  The  anaiytical  potentiai  is  fitted  to  this  data 
base  using  an  iterative  method  invoiving  simultaneous  parameter 
adjustment.  Using  these  methods,  we  have  developed  global 
potential  surfaces  for  the  SiH2^^,  Sl3^^,  S/4^®,  Si2H4^^,  ShHe^^, 
C2H4^o,  and  H2FC-CH2F^^>^^  systems.  In  each  case,  the  global 
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surfaces  correctly  predict  a  wide  array  of  structural,  energetic,  and 
dynamical  data  for  the  system. 


Numerous  homogeneous  gas-phase  reactions  play  a  role  in  silicon 
crystal  growth  via  CVD  methods.  Except  at  very  low  pressures^^,  the  initiating 
reactions  in  the  silicon  CVD  process  involve  the  unimolecular  decomposition  of 
the  starting  material.  If  this  is  silane,  the  principal  reactions  of  interest  are 
(R1)  SiH4  — >  SiHa  +  H  . 

(R2)  SiH4  ™>  SiHa  +  H  ,  followed  by 

(R3)  SiHa  -->  SiH  +  H.  and 

(R4)  SiHa  — >  Si  +  Ha  . 

Once  silicon  atoms  are  formed  via  reaction  R4,  clustering  processes  are 
initiated.  Silicon  dimers  may  be  formed  directly  by  termolecular  reactions  such 
as 

(R5)  Si  +  Si  +  M  — >  Sia  +  M  . 

where  M  is  any  third  body  or  via  polymer  decomposition  pathways  involving 
higher  silanes: 

(R6)  SiH4  +  SiHa  — >  SiaHe  , 

(R7)  SiaHs  — >  SiaH4  +  Ha  , 

(R8)  SiHa  +  SiHa  ™>  SiaH4  ,  followed  by 

(R9)  SiaH4  — >  SiaHa  +  Ha  ,  and 

(R10)  SiaHa  — >  Sia  +  Ha- 

Once  the  Sia  dimer  is  formed,  higher  silicon  clusters  may  be  formed  via  a 
series  of  monomer  or  dimer  addition  reactions  of  the  form 
(R11)  Si  +  Sin  -->  Sin+1  and 

(R12)  Sia  +  Sin  -->  Sln+2  . 

provided  the  adducts,  Sin+i  and  Sin+2  >  can  be  stabilized  by  collisions  with  a 
third  body  before  they  undergo  unimolecular  decomposition.  Coltrin,  Kee,  and 
Miller23  and  Breiland,  Coltrin,  and  Ho^^  have  suggested  that  reactions  R6-R10 
are  the  major  mechanistic  pathways  leading  to  silicon  dimer  formation. 
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The  formation  of  such  clusters  are  a  likely  source  of  particulate  formation 
in  silicon  CVD23, 25-27  Calculations  by  Coitrin  et  al.^^  have  predicted  that  Si2 
and  Sis  make  the  dominant  contribution  to  the  overall  deposition  rate  in  silicon 
CVD  from  silane  at  temperatures  above  1 150°  C  when  the  carrier  gas  is  helium; 
these  results  for  Sis  in  a  silane  CVD  reactor  have  received  qualitative 
experimental  confirmation^^-  Consequently,  our  research  group  has  devoted 
considerable  attention  to  the  study  of  the  rates  and  mechanisms  of  such 
clustering  reactions. 

If  silicon  CVD  is  done  using  disilane  as  the  starting  material,  then  the 
initial  homogeneous  reactions  are  much  more  complex.  Possible  channels  of 
importance  include 

(R13)  SiaHe  --->  SiH4  +  SiHa  , 

(R14)  SiaHs  --->  2SiH3, 

(R15)  SiaHe  — >  SiaHs  +  H  , 

(R16)  SiaHe  — >  H3Si=SiH  -»■  Ha  , 

(R17)  SiaHs  — >  HaSi=SiHa  +  Ha  , 

(R18)  SiaH4  — >  HaSi=SiH  H, 

(R19)  SiaH4  — >  HaSi=Si  +  Ha  , 

(R20)  SiaH4  — >  HSi=SIH  +  Ha  ,  and 

(R21)  SiaH4  — >  2SiHa. 

To  date,  we  have  conducted  theoretical  studies  of  most  of  the  above 
reactions.  These  investigations  have  employed  quasiclassical  trajectory 
techniques^  and  three  classical  variational  phase-space  methods.  The  first  of 
the  variational  phase-space  methods^s  uses  an  RRKM  formalism  but  replaces 
the  quantum-state  counts  with  analogous  classical  sums  over  phase  space. 
The  calculated  rate  is  then  variationally  minimized  with  respect  to  the  location  of 
the  critical  surface.  This  procedure  is  based  upon  the  work  of  Farantos,  Murrell, 
and  Hajduk29,  Bhuiyan  and  Hase^o,  and  that  of  Noid,  Koszykowski,  Tabor,  and 
Marcus^i  who  have  shown  that  there  is  often  excellent  agreement  between  the 
exact  quantal  sums  and  the  classical  phase-space  integrals.  The  second 
variational  method32  (MCVPST)  is  an  application  of  Doirs4.33  generalization  of 
Slater  theory  in  which  he  has  replaced  the  averaging  over  the  normal  mode 
vibrational  phases  with  an  average  over  the  microcanonical  ensemble.  Doll 
used  this  method  to  examine  the  dissociation  of  four-atom  clusters23  and 
Adams34  employed  similar  procedures  to  obtain  vibrational  predissociation 
rates  for  van  der  Waals  clusters.  Finally,  we  have  adapted  the  Efficient 


Microcanonical  Sampling  (EMS)  methods  developed  by  Hippier  et  al.^  and 
Nyman  et  to  Doll's  method^^  and  coupled  this  with  analytical  evaluation  of 
the  integrals  over  the  momentum  spacers  to  produce  a  highly  efficient 
procedure  for  the  study  of  unimolecular  decomposition  reactions  in  very 
complex  molecular  systems.  We  denote  this  procedure  by  EMS-TST. 

We  have  also  extended  the  MCVPST  method  to  permit  calculation  of 
microcanonical  unimolecular  decomposition  rates  out  of  specified  initial 
rotational  states^®.  This  method  is  a  variant  of  that  originally  developed  by 
Wardlaw  and  Marcus®^  who  computed  such  rates  using  a  transition-state 
procedure  in  which  the  "transitional"  modes  are  treated  with  Monte  Carlo 
methods  while  the  quantum  states  for  the  remaining  degrees  of  freedom  are 
counted  directly.  Application  of  the  rotational-state  specific  version  of  the 
MCVPST  method  to  the  decomposition  of  CH4  on  the  ab  initio  potential  surface 
developed  by  Duchovic,  Hase,  and  SchlegeP®  gave  results  in  accord  with 
previous  findings  by  Hase  and  Wolf®®  and  by  Miller  and  Brown^o  and  showed 
that  the  method  was  more  than  an  order  of  magnitude  more  computationally 
efficient  than  the  corresponding  classical  trajectory  calculations. 

The  dynamics  of  the  silicon  clustering  reactions,  R5,  R11,  and  R12, 
leading  to  the  formation  of  dimers,  trimers,  and  tetramers  have  been  examined 
using  trajectory  methodsi^-^®-'^T42.  jhe  potential-energy  hypersurfaces  for  the 
Sia  and  Si4  systems  have  been  obtained  using  the  methods  described  in 
Section  II. B.  The  data  bases  required  for  the  surface  fitting  were  obtained  from 
the  ab  initio  results  reported  by  Peyerimhoff  and  Buenker^®  for  the  Sia  dimer 
and  by  Raghavachari^^-^s  for  the  Sia  and  SU  clusters.  The  Sia  calculations'^® 
were  done  at  the  C!  level  while  Raghavachari'^'^-'^®  carried  out  both  Cl  and  MP4 
calculations  using  (6s,5p.2d,1f)  basis  sets  for  the  trimer  and  MP4  calculations 
with  a  6-31 G*  basis  set  for  the  tetramer.  The  results  provide  data  for  several 
Sia  and  Si4  electronic  states  having  various  geometries.  In  addition,  we  have 
incorporated  the  experimental  data  reported  by  Chatillon  ef  a/.'^®  and  Huber  and 
Herzberg47.  For  the  Sia  systemic,  the  calculated  rms  deviation  of  the  energies 
obtained  from  the  global  surface  from  the  ab  initio  results'*'*  is  0.117  eV. 
Equilibrium  geometries  and  energies  are  essentially  in  exact  accord  with  the 
electronic  structure  calculations.  The  fundamental  normal  mode  frequencies 
predicted  by  the  global  surface  are  within  37  cm-i  or  less  of  those  given  by  the 
MP4  calculations^'*.  Results  of  similar  quality  are  obtained  for  the  tetramer*®. 
In  this  case,  the  average  difference  between  the  predictions  from  the  global 


surface  and  the  ab  initio  data'*^  for  the  internuclear  distances  for  the  minimum 
energy  square,  tetrahedron,  and  linear  Si4  structures  is  0.166  a.u.;  the  error  in 
any  of  these  distances  is  less  than  5%.  The  average  deviation  of  the  calculated 
fundamental  frequencies  for  the  ground-state  Si4  rhombus  structure  is  35.1 
cm'T 

The  rates  and  dynamics  for  reaction  R5  with  M  =  Ar  or  Si  have  been 
investigated  at  temperatures  of  800,  1000,  and  1200  The  calculated 

rate  coefficients  for  the  [Ar,Si,Si]  system  all  lie  in  the  range  of  1.34-1.46  x  10^® 
cm®/mo|2-s.  If  rotationally  trapped  dimers  are  included,  the  results  are  in  the 
range  2.51-2.68  x  10^®  cm^/mol^-s.  The  weak  temperature  dependence  is 
characterized  by  an  activation  energy  of  1 .2  kcal/mol.  When  Si  is  the  third  body, 
the  rates  are  more  than  an  order  of  magnitude  larger.  Four  mechanistic 
pathways  leading  to  recombination  have  been  identified.  These  are  (a)  direct 
energy  exchange,  (b)  direct  atom  exchange,  (c)  complex  formation,  and  (d) 
metastable  formation  due  to  a  rotational  barrier.  For  the  [Si,Si,Si]  system,  the 
relative  contribution  of  these  pathways  are  [(a)+  (b)]  {66%},  (c)  (6%),  and 
(d){28%}.  Internal  energy  distributions  for  product  Sia  dimers  have  also  been 
obtained  for  both  systems, 

Silicon  clustering  via  reaction  reaction  R1 1  with  n  =  2  and  3  have  also 
been  examined^ 8  Formation  cross  sections  for  the  trimer  are  found  to 
exhibit  a  near  double  exponential  dependence  upon  relative  translational 
energy.  Calculated  thermal  rate  coefficients  are  on  the  order  of  10^5  cm^/mol-s 
and  show  a  negative  temperature  dependence  at  temperatures  below  200  K. 
Microcanonical  rate  coefficients  at  six  different  energies  have  been  calculated. 
These  results  follow  an  RRK  energy  dependence  with  s  =  2.67.  The  average 
Sis  lifetimes  are  found  to  lie  in  the  range  0.47  -  50.2  ps.  The  formation  cross 
section  for  Si4  peaks  sharply  near  Et  =  0  and  falls  off  in  near  linear  fashion  with 
energy.  The  thermal  rate  coefficients  for  tetramer  formation  lie  in  the  range  6-8 
X  10^^  cm3/mol-s  for  the  temperature  range  800-1500  K.  Microcanonical  rate 
coefficients  are  well-described  by  an  RRK  expression  with  s  =  4.67.  The 
average  Si4  lifetime  is  computed  to  be  413  ps  at  800  K.  This  lifetime,  as  well  as 
that  for  the  trimer,  are  not  sufficiently  long  to  permit  a  stabilizing  collision  to 
occur  at  pressures  characteristic  of  low-pressure  CVD  experiments. 
Consequently,  we  may  conclude  that  if  Sia  and  Si4  play  a  role  in  silicon  surface 
growth  in  CVD  experiments,  their  formation  must  involve  a  pathway  other  than 
reaction  R1 1.  Some  possible  pathways  are  discussion  in  Section  III.B. 


The  decomposition  pathways  and  rates  for  silane,  reactions  R1-R4,  have 
been  investigated  in  the  initial  phases  of  the  AFOSR  project28.32,48  During  the 
first  two  years  of  the  current  grant,  we  have  directed  our  attention  toward  the 
study  of  the  more  complex  Si2H4  and  Si2H6  systems,  reactions  R13- 
R21 16,19,49. 

The  potential-energy  surfaces  for  Si2H4  and  Si2H6  were  developed 
using  the  methods  of  Section  II.B.  The  ab  initio  calculations  reported  by  Ho, 
Coltrin,  Binkley,  and  Melius^O-Si,  Gordon  and  co-workers52,53  Binkley54,  and 
Peyerimhoff  and  Buenker^  along  with  important  experimental  data  obtained  by 
Berkowitz  et  alP^,  JasinskiS®,  and  Inoue  and  Suzuki57  make  up  the  data  base 
to  which  the  global  surfaces  are  fitted.  These  electronic  structure  calculations 
use  Cl  and  MP4  methods  to  obtain  bond  enthalpies,  equilibrium  geometries, 
fundamental  vibrational  frequencies,  and  transition-state  enthalpies  for  SiHn 
and  Si2Hn  (n  =  0-6)  molecules.  The  experimental  data  provide  Si-H  bond 
energies  and  activation  energies  for  some  of  the  important  decomposition 
channels.  The  resulting  global  surface  for  disilane  contains  40  many-body 
terms  with  88  parameters  that  describe  the  geometries,  energies,  and 
stretching,  bending,  and  torsional  motions  for  Si2H6,  Si2H5  ,  H3Si=SiH, 
H2Si=SiH2,  H2Si=SiH,  H2Si=Si,  HSi=SiH,  H2,  SiH2,  and  Si2  species  as  well 
as  the  barrier  heights  and  reaction  profiles  for  all  of  the  important  reaction 
channels  of  the  system  described  by  reactions  R13-R21.  In  general,  the 
equilibrium  bond  lengths  and  angles  given  by  this  global  surface  are  in 
agreement  with  the  ab  initio  results  to  within  0.03  A  and  0.5°,  respectively.  The 
calculated  heats  of  reaction  for  24  different  channels  are  in  excellent  agreement 
with  the  ab  initio  MP4  calculations50-54  and  with  the  experimental  data55-57 
The  average  error  is  2.79  kcal/mol.  The  average  deviation  of  the  predicted 
fundamental  frequencies  for  Si2H6,  Si2H5,  H3Si=SiH,  H2Si=SiH2,  H2Si=SiH, 
H2Si=Si,  and  SiH2  from  the  results  reported  by  Ho  ef  a/.50.5i  js  55.7  cm-i.  The 
computed  barrier  heights  are  in  accord  with  measured  thermal  activation 

energies56,57 

The  unimolecular  decomposition  reactions  of  Si2H4  and  Si2H6,  reactions 
R13-R21,  have  been  investigated  on  the  above  global  potential-energy 
surfaces  using  trajectory2. 16,19,49  and  variational  EMS-TST  methods25,58.  in 
the  latter  calculations,  the  microcanonical  rate  coefficient,  k(E),  for  unimolecular 
reactions  is  expressed  as  an  average  over  the  microcanonical  ensemble^s 


(2)  k(E)  =  0.5  /  dr5[H(r)-E]  5[qRc  -  qcl  IVrcI  /  /  cir5[H(r)-El , 


where  r  is  the  complete  set  of  position  and  momentum  coordinates  [q,  p],  H(r) 
is  the  Hamiltonian  of  the  system  excluding  center-of-mass  motion,  qRc  =  qRc(q) 
is  the  reaction  coordinate,  which  may  be  a  function  of  some  or  all  of  the 
coordinates  q,  qc  is  the  critical  value  required  for  reaction,  and  |Vrc|  is  the 
absolute  value  of  the  velocity  component  perpendicular  to  the  critical  surface. 
The  integrals  over  r  in  Eq.(2)  are  understood  to  be  over  the  reactant  portion  of 
phase  space.  Equation  (2)  is  equivalent  to  the  microcanonical  rate  coefficient  of 
generalized  transition-state  theory^^  for  the  special  case  that  the  critical  surface 
separating  reactant  and  product  configurations  is  a  function  of  coordinates  only. 

In  the  MCVPST32  application  of  Eq.{2),  the  integrals  are  evaluated  using 
Metropolis  sampling  over  the  entire  6N-dimensional  phase  space  r.  In 
addition,  the  delta  function,  8[H(r)-E],  is  replaced  with  a  prelimit  form.  While 
this  procedure  serves  to  increase  the  rate  of  convergence  of  the  Monte  Carlo 
sums  that  are  used  to  approximate  the  integrals,  it  allows  points  lying  off  the 
energy  shell  to  make  unwanted  contributions  to  the  rate  coefficient. 

In  the  present  calculations^s.ss^  vve  have  introduced  methods  that 
obviate  both  of  the  above  problems.  By  incorporation  of  the  Efficient 
Microcanonical  Sampling  (EMS)  procedures*^,  we  may  eliminate  the  need  to 
employ  a  prelimit  form  of  the  delta  function  and  thereby  remove  any  contribution 
to  k(E)  from  phase-space  points  that  lie  off  the  energy  shell.  In  addition,  by 
exploiting  the  separability  between  potential  and  kinetic  energy  in  the 
Hamiltonian,  wo  may  execute  the  integrations  over  the  momentum  coordinates 
in  Eq.(2)  analytically  and  thereby  derive  a  closed-form  expression  for  the 
average  velocity  perpendicular  to  the  critical  surface  for  the  case  of  simple  bond 
fission  reactions^s.  When  this  is  done,  Eq.(2)  becomes 

(3)  k(E)  =  0.5  <1Vrc|>  I  drW(q)6[qRc  -  qc]  /  1  drW(q)  , 

where  W(q)  is  the  appropriate  EMS  statistical  weight  that  has  been  shown  to 
be®*^ 

(4)  W(q)=  [E-V(q)](3N-5)/2  . 

For  bond  fission  reactions,  we  have  shown^s  <|VRci>  to  be  given  by 


(5)  <|VrcI>  =  [2K/n4R]1/2[(3N-5)/2!]/[(3N-4)/2]!  . 


where  |ir  is  the  reduced  mass  associated  with  the  reaction  coordinate  and 

(6)  K  =  E  -  V(qc) . 

To  increase  the  rate  of  convergence  of  Eq.(3),  we  have  also  included 
importance  sampling  to  increase  the  probability  that  the  Markov  walk  samples 
regions  of  higher  potential's.  To  do  this,  we  have  employed  an  importance 
sampling  weight,  l(q),  with  the  same  form  as  W(q): 

(7)  l(q)  =  [E  -  V(q)]a . 

which  gives  an  effective  weight  function  for  the  Markov  walk  of 

(8)  W(q)/I(q)  =  [E  -  V(q)]{(3N-5)/2}  -  a . 

The  exponent  a  is  adjusted  to  provide  a  satisfactory  Markov  walk  for  accurate 
evaluation  of  both  the  numerator  and  denominator  of  Eq.(3). 

The  recombination  of  silylene  molecules  to  form  disilene, 

(R22)  SiHa  +  SiHa  --->  SiaH4  , 

has  been  examined  using  trajectory  methods^®.  Cross  sections  and  thermal 
rate  coefficients  have  been  determined  over  the  range  300-1500  K.  Thermally 
averaged  formation  cross  sections  vary  from  66.3  to  28.7  A2  over  this  range. 
The  corresponding  thermal  rate  coefficients  lie  in  the  range  2-4  x  10"''^  cm^/mol- 
s.  We  have  found  that  the  reaction  exothermicity  is  primarily  partitioned  into  the 
Si-Si  stretch  and  the  H-Si-H  bending  modes  upon  formation  of  SiaH4. 

To  date,  there  have  been  no  measurements  of  the  rate  coefficients  or 
cross  sections  for  reaction  R22.  However,  some  photolysis  data  have  been 
reported  for  the  analogous  reaction  of  SiHa^O’^^, 

(R23)  SiHa  +  SiHa  ™>  SiaHs  • 

The  reported  rate  coefficient  for  R23  was  0.5  x  10^^  cm^/mol-s  at  373  K.  This 
value  is  1/4  the  corresponding  rate  computed  for  reaction  R22.  Since  the 
attractive  interaction  and  Si-Si  bond  energy  is  greater  for  the  recombination  of 
SiHa  than  for  SiHa,  and  since  we  would  expect  more  steric  hindrance  in  the 
recombination  process  for  SiHa,  our  computed  results  seem  to  be  in  reasonable 
accord  with  the  experimental  data. 

The  dissociation  dynamics  of  SiaH4  have  been  examined  in  great  detail 
over  the  internal  energy  range  5.0  <  E  <  9.0  eV^9.  In  this  range,  the  most 
important  dissociation  channel  is  three-center  elimination  of  molecular 


hydrogen  leading  to  H2Si=Si,  reaction  R19.  For  energies  below  7.0  eV,  the 
other  decomposition  channels  are,  in  order  of  importance,  Si-Si  bond  rupture 
(R21),  four-center  Ha  elimination  (R20),  and  Si-H  bond  rupture  (R18).  At  or 
above  8.0  eV,  four-center  Ha  elimination  replaces  Si-Si  bond  rupture  as  the 
second  most  important  dissociation  channel.  Microcanonical  rate  coefficients 
for  all  of  these  channels  have  been  obtained.  We  find  that  three-center  Ha 
elimination  involves  a  simultaneous  rupture  of  both  Si-H  bonds  whereas  the 
four-center  elimination  channel  is  found  to  proceed  by  a  hydrogen-atom  transfer 
process  followed  by  Ha  elimination.  Except  for  a  propensity  to  form  Ha  with 
excess  rotational  energy,  the  energy  partitioning  among  the  products  is  nearly 
statistical. 

The  SiaHe  dissociation  dynamics  at  seven  internal  energies  ranging  from 

5.31  to  9.31  eV  have  been  determined  using  both  trajectory  and  EMS-TST 
methods^9'35.58.  At  low  internal  energies,  dissociation  leading  to  SiH4  +  SiHa, 
reaction  R13,  dominates  the  dynamics.  For  energies  in  the  range  5.31  <  E  < 

6.31  eV,  the  various  decomposition  channels  are,  in  order  of  importance, 
hydrogen-atom  transfer  leading  to  SiH4  +  SiHa  (R13),  Si-Si  bond  rupture  (R14), 
three-center  Ha  elimination  (R16),  Si-H  bond  rupture  (R15),  and  four-center  Ha 
elimination  (R17).  At  higher  energies,  entropy  effects  cause  an  inversion  of  this 
ordering  such  that  Si-Si  bond  rupture  becomes  the  major  reaction  channel 
followed  by  three-center  Ha  elimination,  SiH4  +  SiHa  formation,  Si-H  bond 
rupture,  and  four-center  Ha  elimination.  Our  results  indicate  that  formation  of 
disilene  in  silicon  CVD  experiments  occurs  predominately  via  the  formation  of 
HaSisaSiH  from  three-center  Ha  elimination  followed  by  a  low-barrier  hydrogen- 
atom  transfer  process.  For  all  decomposition  channels,  most  of  the  available 
energy  is  partitioned  into  vibrational  modes  of  the  products.  To  a  large  extent, 
product  energy  partitioning  is  found  to  be  governed  by  statistical  considerations. 
We  find  that  while  three-center  Ha  elimination  is  a  concerted  reaction,  it 
probably  does  not  proceed  via  a  symmetric  pathway.  Reaction  R13  is  found  to 
concerted,  but  four-center  Ha  elimination  involves  the  rupture  of  one  Si-H  bond 
followed  by  a  hydrogen  transfer  and  a  subsequent  Ha  abstraction  reaction  to 
give  the  Ha  +  SiaH4  products. 

One  of  the  most  interesting,  and  perhaps  important,  findings  of  these 
studies  is  the  fact  that  the  EMS-TST  calculations^s  of  the  microcanonical  rate 
coefficients  for  Si-H  and  Si-Si  bond  rupture  in  SiaHe  demonstrate  the  existence 
of  pronounced  non-statistical  effects  in  spite  of  the  fact  that  the  total  internal 


energy  content  is  well  in  excess  of  the  dissociation  thresholds.  We  have  found 
evidence  that  such  non-statistical  effects  may  be  widespread  among  polyatomic 
molecules.  These  results  call  into  question  the  practice  of  using  statistical 
methods,  such  as  transition  state  theory,  to  investigate  these  reactions. 

Summary: 

We  have  employed  trajectory  and  variational  phase-space 
methods  to  examine  the  dynamics  of  silicon  clustering  reactions, 
silylene  recombination,  and  unimolecular  dissociation  reactions  of 
SIH4,  SiH2,  Si2H4,  and  Si2He.  Global  potential-energy  surfaces 
are  developed  for  each  of  these  systems  using  the  results  of 
electronic  structure  calculations  and  experimental  measurements 
as  the  data  base  to  which  we  fit  many-body  functional  forms 
suggested  by  physical  and  chemical  considerations.  The  global 
potentials  are  shown  to  give  good  to  excellent  results  when 
compared  to  measured  and/or  calculated  equilibrium  geometries, 
energies,  fundamental  frequencies,  and  barrier  heights.  Dynamical 
results  Include  cross  sections,  thermal  and  microcanonical  rate 
coefficients,  lifetimes,  detailed  mechanisms,  and  energy 
partitioning  information  for  all  open  reaction  channels.  Striking 
non-statistical  effects  are  found  for  the  bond  fission  reactions  of 
disllane. 

IIC.2.  Heterogeneous  Reactions 

The  surface  processes  of  interest  in  silicon  CVD  include  inelastic  surface 
scattering  and  sticking,  dissociative  chemisorption,  surface  diffusion,  tunneling, 
and  chemical  reactions.  To  investigate  these  processes,  we  have  employed 
MCVPST32,  trajectory®^,  and  master  equation®^  methods  to  examine  surface 
diffusion,  and  WKB  semiclassical  methods  coupled  with  MCVPST  calculations 
to  study  the  effects  of  tunneling  on  hydrogen-atom  diffusionf^.  inelastic 
scattering  events  and  sticking  probabilities  have  been  computed  using 
semiclassical  wavepacket®®,  perturbation®-®,  and  trajectory®®  methods.  To 
date,  the  surface  reactions  investigated  include  the  dissociative  chemisorption 
of  H2  on  Si(111)  surfaces  with  and  without  partial  hydrogen  coverage®^-®®, 
thermal  desorption  of  molecular  hydrogen  from  Si(111)  surfaces®®-^®,  the 


chemisorption  and  subsequent  decomposition  of  SiHa  on  Si(111)^^  and 
reconstructed  Si(111)-(7x7)  surfaces^^^  and  the  chemisorption  of  silicon 
clusters,  Sin  (n  =  2,3,4),  on  Si(1 1 1)-(7x7)  surfaces^3.  Of  particular  interest  in 
these  studies  are  the  effects  of  surface  temperature,  surface  orientation  and 
reconstructions,  surface  dislocations  and  kinks,  and  surface  coverage. 

The  problem  of  obtaining  sufficiently  accurate  potential-energy  surfaces 
is  even  more  acute  for  heterogeneous  systems  than  is  the  case  for  the  gas- 
phase  reactions  discussed  above.  Except  for  a  few  cluster  calculations  for 
selected  systems,  along  with  some  thermodynamic  and  spectroscopic  data, 
there  is  usually  very  little  information  available  on  the  topographical  features  of 
these  potential  surfaces.  Consequently,  we  must  content  ourselves  with  the  use 
of  more  empirical  surfaces.  Although  we  might  expect  such  hypersurfaces  to 
yield  at  most  qualitative,  or  in  some  cases,  semi-quantitative  information,  we 
have  found  that  the  results  can  be  surprisingly  good  provided  care  is  taken  to  fit 
the  known  energetics  of  the  system  as  well  as  possible.  Several  examples 
taken  from  our  recent  studies  are  given  below.  In  spite  of  some  notable 
successes  to  the  contrary,  we  generally  regard  our  theoretical  results  in 
heterogeneous  systems  to  be  only  indicators  of  the  true  system  behavior,  which 
is  often  unknown.  Such  "indicators"  are  useful  in  themselves  and  often  point 
the  way  to  new  experiments  or  theoretical  methods  that  will  provide  more 
accurate  information. 

In  our  initial  studies,  we  have  investigated  inelastic  scattering,  sticking, 
and/or  dissociative  chemisorption  of  Si,  H,  H2,  and  SiH2  incident  on 
unreconstructed  Si(100)  and  Si(111)  surfaces  and  partially-to-fully  covered 
(with  H  atoms)  Si(111)  surfaces®5-68,7i,72  7^0  lattice  models  used  in  these 

calculations  have  9  to  60  atoms  in  the  primary  reaction  zone  (P-zone)  and  70  or 
more  in  the  secondary  Q-zone^'*.  Both  Weber’s  parametrization^s  of 
Keating's^s  potential  and  a  dissociative  potential  developed  by  Brenner  and 
Garrison77  were  used  to  represent  the  lattice  interactions.  To  incorporate  the 
effects  of  relaxation  of  the  P-zone  atoms  to  the  bulk  phonon  modes,  we 
developed  a  velocity  reset  method®®  using  concepts  analogous  to  those 
employed  by  Riley,  Coltrin,  and  Diestler^®.  In  this  method®®,  the  velocities  of 
the  peripheral  atoms  of  the  P-zone  are  reset  at  periodic  intervals  according  to  a 
prescription  based  on  the  local  temperature  and  the  temperature  of  the  bulk. 
Berendsen  et  al7^  have  shown  that  this  method  is  equivalent  to  a  Langevin 
procedure  under  conditions  that  all  frictional  coefficients  are  equal  and  there  is 


no  Gaussian  white  noise  term.  In  our  more  recent  studies,  we  have  utilized  the 
Riley  et  al.^^  procedure  to  incorporate  bulk  relaxation  effects. 

The  principal  results  obtained  from  the  above  studies  are  as  follows:  (1) 
Sticking  coefficients  for  Si  atoms  and  SiHa  molecules  are  near  unity  at  all 
surface  temperatures.  (2)  The  major  processes  observed  in  collisions  of 
hydrogen  atoms  with  partially  covered  Si(111)  surfaces  are  direct  scattering, 
exchange  scattering,  direct  chemisorption,  and  chemisorption  with  exchange. 
(3)  For  the  case  of  H-atom  adsorption  on  partially-covered  Si(1 1 1 )  surfaces,  the 
extent  of  surface  coverage  is  found  to  have  only  a  small  effect  upon  the 
calculated  sticking  probabilities.  (4)  Energy  transfer  to  the  bulk  phonon  modes 
is  found  to  be  a  first-order  process  in  which  energy  is  transferred  mainly  from 
stretching  motions  of  the  adsorbate  to  the  bending  modes  of  the  silicon  surface. 
Rate  coefficients  for  this  transfer  decrease  with  increasing  surface  temperature 
and  are  on  the  order  of  1-2  x  10^2  5-1  .  (5)  The  nature  of  H2  dissociative 
chemisorption  on  Si(111)  surfaces  is  qualitatively  similar  to  that  observed  by 
Wolken  and  co-workers^3-82  for  the  H2-W(100)  system  and  by  Lee  and 
DePristo®3.84  for  H2-Ni  systems.  The  quantitative  differences  are  mainly  due  to 
the  increased  exothermicity  of  H2  chemisorption  on  Si(111)  surfaces.  H2 
adsorption  in  this  system  is  always  accompanied  by  dissociation  and  the 
sticking  probabilities  are  virtually  independent  of  surface  relaxation,  the  size  of 
the  lattice  model,  and  the  nature  of  the  lattice  potential.  Si-H  vibrational  energy 
transfer  rates  to  the  lattice  are  increased  and  surface  mobilities  are  decreased 
by  the  inclusion  of  relaxation  to  bulk  modes  but  the  magnitude  of  these  effects  is 
less  than  a  factor  of  two. 

During  the  grant  period  for  AFOSR-89-0085,  we  have  directed  our 
attention  toward  studies  of  the  chemisorption  and  subsequent  reaction  of 
SiH2^2^  Si2,  Si3,  and  Si4  clusters,  and  silicon-atom  diffusion  rates®5.86  on 
reconstructed  Si(1  ll)-(7x7)  surfaces.  We  have  also  developed  perturbation 
methods  that  significantly  reduce  the  computational  requirements  for  the  study 
of  surface  processes®-®. 

The  dynamics  of  chemisorption  and  decomposition  of  SiH2  on  Si(111)- 
(7x7)  surfaces  have  been  investigated  using  trajectory  methods72.  The  Binnig 
et  al.^^  model  is  employed  to  represent  the  (7x7)  reconstruction.  The  lattice 
model  in  these  studies  consists  of  the  12  silicon  adatoms  found  by  Binnig  et 
al.^’^  in  the  topmost  layer  of  the  surface,  81  silicon  atoms  in  layer  B,  99  in  layer 
C,  and  99  in  layer  0  for  a  total  of  291  atoms.  In  some  calculations,  we  have 


employed  the  Keating  lattice  potential^®  with  the  force  constants  modified  to 
yield  the  lattice  strain  energy  predicted  by  the  Brenner-Garrison  potential^^.  In 
other  studies,  this  lattice  potential  has  been  replaced  with  one  suggested  by 
Bolding  and  Andersen®®  to  describe  silicon  clusters,  crystals,  and  surfaces. 
This  potential  is  a  generalization  of  the  bonding  formulation  used  by  Tersoff®®  in 
which  the  total  potential  is  represented  by  a  sum  of  two-body  terms  each  of 
which  contains  an  environment-dependent  "interference"  function.  The  SiHa- 
surface  interaction  potential  is  fitted  to  calculated  bond  energies,  structures, 
measured  and/or  calculated  vibrational  frequencies,  and  to  ab  initio 
calculations  of  the  repulsive  interactions  between  SiHe  and  the  closed-shell 
lattice  atoms  using  methods  similar  to  those  described  in  Section  II. B. 

Our  results  predict  SiH2  sticking  probabilities  to  be  near  unity  on  the 
reconstructed  surface.  Subsequent  to  chemisorption  the  major  mode  of  SiHa 
decomposition  is  found  to  be  direct  release  of  molecular  hydrogen  into  the  gas 
phase; 

(R23)  *-SiH2  — >  *-Si  +  H2  . 

where  *  represents  a  surface  binding  site.  Hydrogen-atom  dissociation  to 
adjacent  lattice  sites  is  found  to  be  the  secondary  decomposition  channel: 

(R24)  *-SiH2  +  *  ->  *-SiH  +  *-H  . 

The  calculated  rate  coefficients  are  3.4  x  10^®  and  0.8  x  10^®  s'''  for  reactions 
R23  and  R24,  respectively. 

At  the  time  of  publication  of  the  above  results,  the  only  experimental  data 
related  to  the  surface  decomposition  of  SiH2  was  that  obtained  by  Farnaam  and 
Olander®®  from  their  molecular  beam  studies  of  silane  decomposition  on 
Si(1 1 1)-(7x7)  surfaces.  They  proposed  that  SiH4  chemisorption  occurs  via  the 
reaction 

(R25)  SiH4  +  *  +  *-*  — >  *-SiH2  +  H-*-H  +  *  . 

Subsequent  to  chemisorption,  these  investigators  observed  only  a  single 
decomposition  channel  for  chemisorbed  SiH2,,  the  concerted  release  of 
molecular  hydrogen,  reaction  R23.  The  measured  thermal  rate  coefficient  for 
R23  at  1000  K  is  5.9  x  10^  s'^  or  less®®.  The  fact  that  experiments  involving 
simultaneous  surface  reaction  of  SiH4  and  S1D4  led  only  to  the  formation  of  H2 
and  D2  with  no  HD  being  produced,  led  Farnaam  and  Olander®®  to  suggest  that 
decomposition  via  hydrogen-atom  transfer  to  adjacent  lattice  sites,  reaction 
R24,  does  not  occur. 


In  order  to  compare  our  computed  rates  for  reaction  R23  with  the 
molecular  beam  data^o,  we  must  extrapolate  to  the  thermal  regime.  In  the 
Farnaam-Olander  experiments,  *-SiH2  is  formed  via  reaction  R25  which  is 
probably  nearly  thermochemically  neutral.  Consequently,  the  measured  rates 
for  H2  release  will  correspond  to  thermal  rates.  In  contrast,  in  our  calculations, 
*-SiH2  is  formed  via  chemisorption  of  SiH2(g)  which  is  highly  exothermic.  As  a 
result,  the  computed  decomposition  rates  correspond  to  that  for  highly  excited 
*-SiH2.  If  we  assume  that  these  rate  coefficients  represent  the  high-temperature 
limiting  values,  we  may  estimate  the  corresponding  thermal  rates  using  the 
calculated  barrier  height  for  reaction  R4  for  the  activation  energy.  The  result  of 
this  calculation  is  a  thermal  rate  coefficient  of  9.4  x  10^  S'1  at  1000  K,  which  is  in 
satisfactory  agreement  with  the  Farnaam-OlanderSO  result  considering  the 
uncertainties  inherent  in  the  required  extrapolation. 

The  major  source  of  disagreement  between  our  model  studies  and  the 
molecular  beam  data^O  involves  the  *-SiH2  decomposition  mechanism. 
Farnaam  and  Olander  find  no  evidence  for  the  occurrence  of  reaction  R24 
whereas  our  calculations  predict  that  it  should  be  a  secondary,  but  easily 
observable,  channel.  These  predictions  have  recently  been  confirmed  by  Gates 
et  a/.9i  using  static  secondary  ion  mass  spectrometry  to  observe  the  silicon 
hydride  species  formed  by  silane  adsorption  on  clean,  single  Si(1 1 1)-(7x7) 
crystals.  Their  data  show  that  *-SiH2  is  the  only  species  formed  on  the  surface 
by  silane  deposition.  Subsequent  to  chemisorption,  the  mass  spectral  peaks 
show  the  presence  of  SiH2  and  SiH  on  the  surface.  Consequently,  reaction 
R24  does  occur  as  predicted  by  our  calculations. 

Most  recently,  we  have  initiated  studies  of  the  chemisorption  and  surface 
reaction  dynamics  of  silicon  clusters.  Sin  (n  =  2,3,4),  on  Si(1 1 1  )-(7x7) 
surfaces^3.  These  model  studies  utilize  the  Binnig  et  surface  described 
above  with  the  Bolding-Anderson^s  lattice  potential.  The  velocity  reset  method 
developed  by  Riley  et  a/.^3  is  employed  to  model  relaxation  to  bulk  phonon 
modes.  To  date,  we  have  calculated  sticking  probabilities  for  each  of  the 
clusters  and  determined  the  various  modes  of  chemisorption.  For  example. 
Figs.  1  and  2  illustrate  the  two  most  important  species  formed  upon 
chemisorption  of  Si4  clusters.  In  the  gas  phase,  Si4  is  a  planar  rhombus. 
Figure  1  shows  that  many  chemisorption  events  lead  to  the  Si4  cluster  being 
bound  by  a  single  Si-Si  bond  with  the  rhomboidal  plane  of  the  cluster  intact  and 
nearly  perpendicular  to  the  surface  plane.  A  second  chemisorption  mode  is 
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shown  in  Fig.  2.  Here,  the  symmetry  of  the  cluster  is  distorted  sufficiently  to 
permit  the  formation  of  two  Si-Si  bonds. 

The  diffusion  rates  of  silicon  atoms  on  Si(111)-{7x7)  surfaces  are  crucial 
factors  in  determining  the  nature  and  rate  of  silicon  crystal  growth.  We  have 
therefore  computed  such  diffusion  rates  on  both  the  Binnig  et  al.^^  and  the 
dimer-adatom-stacking  fault  (DAS)  model  introduced  by  Takayanagai  et 
The  jump  frequencies  for  site-to-site  diffusion  are  computed  using  a  canonical 
MCVPST32  method.  That  is,  Eq.(2)  is  cast  into  the  form 

(9)  F(T)  -gJPoutpv-dS/ J  pdr, 

where  F(T)  is  the  classical  flux  across  the  critical  surface  at  temperature  T,  p  is 
the  phase-space  density  of  points,  S  is  the  critical  surface,  v  is  a  generalized 
velocity  vector,  g  is  a  correction  factor  for  surface  recrossing,  and 

(10)  Pout  =  +1  lor  v-dS  >  0  and  Pout  =  0  otherwise. 

For  a  system  in  thermal  equilibrium,  we  have 

(11)  p  =  po  exp[-BH]. 

Finally,  in  order  to  increase  the  convergence  rate  in  the  Monte  Carlo  evaluation 
of  the  phase-space  integrals  in  Eq.(9),  we  incorporate  importance  sampling. 
Let  Va  be  the  interaction  between  the  lattice  and  the  Si  adatom.  We  define 

(12)  H'  =  H  -  Va. 

Substitution  of  Eqs.(1 1)  and  (12)  into  (9)  yields 

(13)  F(T)  =  g  j  Pout  exp[-BVa]  exp[-6H']  v-dS  /  J  exp[-(3Va]  exp[-BHT  dr  . 

In  practice,  the  integrals  in  Eq.(13)  are  evaluated  using  a  Markov  walk  governed 
by  the  distribution  function  exp[-BH']. 

The  flux  F(T)  will  be  an  upper  limit  to  the  true  classical  rate  as  long  as  the 
system  is  statistical.  We  have  therefore  obtained  close  upper  bounds  for  the 
jump  frequencies  by  minimization  of  F(T)  with  respect  to  the  parameters  of  the 
critical  surface.  The  correction  factor  for  recrossing,  g,  is  computed  by  the 
calculation  of  trajectories  starting  at  the  critical  surface^^.  Once  the  site-to-site 


jump  frequencies  are  in-hand,  diffusion  coefficients  are  obtained  by  solution  of 
the  appropriate  master  equation  to  obtain  the  average  square  displacements  as 
a  function  of  time  which  are,  in  turn,  linearly  related  to  the  diffusion  coefficient®^. 

The  calculated  silicon-atom  diffusion  coefficients  on  the  Binnig  et 
surface  are  D  =  2.15  x  IQ-®  exp[-1.51  eV/kT]  cm2/s.  The  calculated  activation 
energy  of  1 .51  eV  is  in  excellent  accord  with  the  results  obtained  by  Farrow^® 
from  molecular-beam  pyrolysis  data  on  SiH4  deposition.  In  contrast,  diffusion 
on  the  DAS  surface®^  is  predicted  to  be  much  slower  due  to  the  larger  spacing 
between  dangling  bonds  on  this  model  of  the  (7x7)  reconstructions**.  The 
calculated  ratio  D[Binnig  et  a/.]/D[DAS]  is  32  at  1 000  K  and  2.5  x  1 0^  at  300  K. 
In  addition,  we  find  that  preferential  directions  of  diffusion  exist  on  the  Binnig  et 
al.  surface.  These  directions  correspond  to  "gateways"  at  three  of  the  four 
corners  of  the  (7x7)  unit  cell.  In  contrast,  diffusion  on  the  DAS  surface  is 
predicted  to  be  isotropic.  Consequently,  the  extreme  differences  in  rates  and 
directions  of  flow  for  diffusion  on  the  two  models  of  the  reconstruction  suggest 
that  careful  measurement  of  such  quantities  may  provide  an  additional 
experimental  method  for  differentiating  the  two  proposed  models  of  the  surface. 

To  conclude  this  section  on  heterogeneous  processes,  we  briefly 
describe  a  new  perturbation  method  that  we  have  developed  during  the  current 
grant  period  for  the  simplified  treatment  of  gas-surface  interactions®-^.  Our 
method  utilizes  the  usual  separation  of  the  system  Hamiltonian  into  zones.  The 
primary  (P)  zone  generally  includes  those  atoms  directly  involved  in  the  process 
under  consideration  plus  an  arbitrary  number  of  atoms  occupying  nearby  sites. 
The  secondary  or  Q-zone  usually  includes  a  larger  number  of  atoms  further 
removed  from  the  site  at  which  the  process  of  interest  is  occurring.  These  atoms 
are  assumed  to  affect  the  process  only  indirectly  through  the  transfer  of  energy 
in  and  out  of  the  P-zone.  In  some  models,  a  stationary  boundary  zone  (B-zone) 
is  also  included  to  reduce  edge  effects.  The  system  Hamiltonian  is  then  written 

(14)  H  =  Tp  +  Tq  +  Vp  +  VpQ  +  Vq  +  VpB , 

where  Tp  and  Tq  are  the  kinetic  energy  terms  for  the  P-  and  Q-zones, 
respectively.  Vp  and  Vq  contain  the  potential  for  the  P  and  Q-zone  atoms 
alone.  The  interactions  between  P-zone  and  Q-zone  atoms  and  between  P- 
zone  and  boundary  atoms  are  represented  by  Vpq  and  Vpb,  respectively. 


The  various  approximate  methods  that  have  been  suggested  to  date 
differ  primarily  in  the  procedures  used  to  treat  the  above  Hamiltonian.  For 
example,  Riley,  Coltrin,  and  Diestler^^  replace  the  Tq,  Vq,  Vrq,  and  Vrb  terms 
with  a  prescription  that  resets  the  velocities  of  the  P-zone  atoms  that  border  the 
Q-zone  at  specified  time  intervals.  The  most  commonly  used  method  involves 
the  use  of  a  generalized  Langevin  equation  in  place  of  the  Q-zone  terms^Ti2 
In  this  procedure,  the  Tq,  Vq,  Vrq,  and  Vrb  terms  in  Eq.(14)  are  all  replaced  by 
a  set  of  Langevin-type  equations  that  operate  on  the  peripheral  atoms  of  the  P- 
zone.  The  net  result  of  this  replacement  is  an  energy  fluctuation  of  the  P-zone 
atoms  that  represents  the  effect  of  the  "heat  bath"  or  Q-zone  atoms. 

In  the  perturbation  method  we  have  developed,  the  effects  of  the  Q-zone 
atoms  are  incorporated  using  the  unforced  solutions  for  these  atoms.  The 
method  assumes  that  the  Q-zone  motion  is  relatively  unperturbed  by  collision 
events  in  the  P-zone.  The  correct  form  of  the  P-zone-Q-zone  interaction  Vrq  is 
retained  and  the  motion  of  all  P-zone  atoms  is  treated  exactly. 

Let  the  N-atom  system  be  partitioned  such  that  there  are  Nr,  Nq,  and  Nb 
atoms  in  the  P,  Q,  and  B  zones,  respectively.  In  general,  we  will  have  Nr  =  Na  + 
M,  where  Na  is  the  number  of  atoms  in  the  colliding  molecule.  A,  and  M  is  an 

arbitrarily  chosen  number  of  atoms  in  the  lattice.  Let  Uj  (i  =  1,2 . 3Na),  xj  (i  = 

1,2,...,3M),  yi  (i  =  1,2,...,3Nq),  and  Zi  (i  =  1,2,...,3Nb)  denote  the  atomic 
coordinates  for  the  atoms  in  the  incident  molecule,  the  P-zone  lattice  atoms,  and 
the  Q  and  B  zones,  respectively.  In  this  notation,  we  may  write  Eq.(14)  as 

(15)  H  =  TA(Ui)  +  TM(xi)  +  TQ(yi)  +  VRA{uj,Xj)  +  Vr^xj)  +  VQ(yj)  + 

VRQ4xi,yj)  +  VRQA(ui,yi)  +  VRB(uj,Xi:zi)  , 

where  the  Vrb  term  depends  only  parametrically  upon  the  Zj.  In  Eq.(15),  the  P- 
Q  interaction  has  been  written  as  the  sum  of  a  lattice-lattice  interaction, 
VRQ'-(xj,yi),  and  an  interaction  between  the  incident  molecule  and  the  Q-zone, 
VRQ^(ui,yi).  The  P-zone  potential  is  likewise  separated  into  the  sum  of  a  lattice- 
lattice  interaction,  VRL(xi),  and  a  lattice-A  interaction,  VRA(ui,Xi). 

It  is  now  explicitly  assumed  that  the  lattice  Q-zone  motion  is  unperturbed 
by  the  collision  of  molecule  A  with  the  surface.  That  is,  we  assume  that  the 
effect  of  VRQA(ui,yi)  on  the  yi  is  negligibly  small  and  that  the  effect  of  the  lattice- 
lattice  interaction  term,  VRQi-(xj,yj),  upon  the  Q-zone  motion  is  unaffected  by  the 
collision  of  A  with  the  surface.  Linder  these  conditions,  the  Q-zone  motion  may 


be  obtained  by  direct  solution  of  Hamilton's  equations  using  the  lattice 
Hamiltonian 

(16)  HL  =  TM(Xj)  +  TQ(yi)  +  Vpl-(Xi)  +  VQ(yi)  +  VpQ>-(Xj,yi)  . 

The  results  of  such  a  solution  yield  the  Q-zone  atomic  coordinates  as  functions 
of  time.  Direct  substitution  of  these  functions  into  Eq.(15)  gives  a  P-zone 
Hamiltonian  of  the  form 

(17)  Hp  =  TA(Ui)  +  TM(Xi)  +  VpA(uj,X|)  +  VpL(xi)  +  VpQ4xi.yi(t)]  + 

VpQ4ui.yi(t)]+  VpB(Ui,Xi;zi)  . 

It  is  assumed  that  the  P-zone  dynamics  may  be  obtained  by  solution  of 
Hamilton's  equations  using  Hp  as  the  Hamiltonian.  If  the  collision  event  exerts 
only  a  negligibly  small  perturbation  upon  the  Q-zone  motion,  the  results 
obtained  using  Eq.(17)  will  approach  those  resulting  from  a  full  solution  of  the 
complete  problem  using  Eq.(15). 

We  have  applied  the  above  method  to  the  inelastic  scattering  of  NO  from 
Ag(1 11)  surfaces  and  to  chemiso'ption  and  surface  reaction  of  SiHa  on  Si(1 11) 
surfaces8-9.  For  the  NO  system,  the  translational  and  rotational  energy 
distributions  of  scattered  NO  molecules,  the  differential  cross  sections,  and  the 
NO  sticking  probabilities  as  a  function  of  incident  translational  energy  obtained 
from  the  perturbation  calculations  with  Np  =  2  are  in  almost  exact  agreement 
with  results  obtained  from  a  full  treatment  of  the  dynamics.  That  is,  the 
perturbation  method  permits  a  near  exact  determination  of  the  system  dynamics 
from  the  results  of  two-body  trajectories.  Equally  good  results  are  obtained  for 
the  SiH2/Si(1 1 1 )  system  if  the  size  of  the  P-zone  is  increased  to  Np  =  9. 

In  general,  the  accuracy  of  the  perturbation  method  increases  as  the 
incident-to-lattice  atom  mass  ratio  decreases.  A  decrease  in  the  strength  of  the 
interaction  between  the  incident  molecule  and  the  Q-zone,  the  incident 
translational  energy,  or  the  lattice  temperature  also  improves  the  accuracy  of 
the  perturbation  treatment.  The  method  is  therefore  best  suited  to  the  study  of 
inelastic,  light-molecule  collisions  with  heavy-atom  surfaces  at  low 
temperatures.  For  such  cases,  near  exact  results  can  be  obtained  with  very 
small  values  of  Np.  For  strong  interactions  or  for  interactions  of  long  duration, 
the  method  will,  of  course,  fai|9. 


Summary: 

During  the  first  two  years  of  the  current  grant  period,  we  have 
determined  the  reaction  dynamics  of  SiH2  chemisorption  and 
subsequent  surface  reaction  on  reconstructed  Si(111)-(7x7) 
surfaces^^.  The  scattering  and  chemisorption  of  Si2,  Sij,  and  SU 
ciusters^^  and  the  rates  and  patterns  of  silicon-atom  surface 
diffusion  on  these  reconstructed  surfaces  have  also  been 
investigated  using  both  the  BInnig  et  al.^'7  and  the  models 

of  the  (7x7)  surface^^’^^.  The  SiH2  surface  decomposition 
mechanism  has  been  found  to  Involve  two  major  channels,  (a) 
direct  molecular  elimination  of  H2  and  (b)  the  dissociation  of 
hydrogen  atoms  to  adjacent  lattice  sites  to  yield  -SiH  and  -H 
surface  species.  Rate  coefficients  computed  for  each  of  these 
channels  shows  the  first  to  be  the  more  important.  These 
predictions,  although  at  odds  with  older  experimental  data^^,  have 
recently  been  confirmed  by  new  SIMS  measurements^ K  Sticking 
coefficients  and  chemisorption  modes  for  each  of  the  silicon 
clusters  have  been  determined.  Diffusion  rates  and  activation 
energies  for  silicon-atom  diffusion  on  Si(111)-(7x7)  have  been 
obtained.  Diffusion  is  predicted  to  be  non-isotropic  on  the  Binnig 
et  afl^  model  of  the  (7x7)  surface  but  isotropic  on  the  DAS^^ 
surface.  In  addition,  we  predict  very  different  diffusion  rates  on  the 
two  surfaces.  It  is  suggested  that  these  differences  may  provide  an 
additional  means  of  experimentally  differentiating  between  the 
models.  Finally,  we  have  developed  a  new  perturbation  method  for 
treating  gas-surface  interactions  that  significantly  reduces  the 
computational  effort  required^’^.  The  method  is  shown  to  give 
nearly  exact  results  in  favorable  cases. 

II.P.  C^t4_-tLF2  and  the  1.2-Difluoroethane  Systems: 

We  have  carried  out  several  investigations  involving  the  reactions  of  F2 
with  ethylene  and  the  unimolecular  decomposition  dynamics  of  1,2- 
difluoroethane.  The  motivation  for  these  studies  arises  from  two  considerations. 


First,  there  are  very  few  experimental  studies  that  provide  detailed  dynamical 
information  about  the  gas-phase  silicon  reactions.  Consequently,  it  is  difficult  to 
assess  the  accuracy  of  our  theoretical  results  for  these  systems.  In  contrast, 
there  exists  a  wealth  of  experimental  data  related  to  the  above  systems  that 
provide  excellent  benchmarks  that  may  be  used  to  assess  the  accuracy  of  our 
methods.  Secondly,  in  Section  III.D.,  we  propose  to  conduct  a  study  of  chemical 
processes  that  occur  under  conditions  of  close  confinement.  We  intend  to  use 
fluorine  reactions  in  cryogenic  matrices  as  prototypes  of  such  reactions  since 
there  exists  a  large  body  of  experimental  data  on  such  systems.  The  gas-phase 
calculations  will  therefore  serve  as  a  basis  for  the  determination  of  the  effects  of 
close  confinement  upon  these  systems. 

The  principal  reactions  of  interest  for  the  ethylene  +  F2  system  are 


(R26) 

C2H4  +  F2  — > 

[H2FC-CFH2]*  • 

(R27) 

C2H4  +  F2  -> 

H2C-CH2F  +  F ,  and 

(R28) 

C2H4  +  F2  — > 

H2C=CHF  +  HF. 

Once  1 ,2-difluoroethane  is  formed,  the  mechanism  and  associated  rates  for  its 
decomposition  become  the  focus  of  attention.  This  decomposition  reaction  may 
proceed  via  several  channels  that  include 


(R29) 

H2FC-CH2F  -> 

H2C=CHF  +  HF 

(R30) 

H2FC-CH2F  -> 

H2FC-CH2  +  F 

(R31) 

H2FC-CH2F  — > 

H2FC-CHF  +  H 

(R32) 

H2FC-CH2F  -> 

2  CH2F  , 

A  global  potential-energy  surface  for  the  C2H4F2  system  has  been 
developed  using  the  methods  described  in  Section  II.B.20.21  The  predicted 
barrier  heights  for  four-center  F2  addition  across  the  ethylene  double  bond 
(R26),  for  HF  elimination  from  1 ,2-difluoroethane  (R29),  for  hydrogen-atom 
abstraction  by  a  fluorine  atom,  and  for  rotation  about  the  C=C  double  bond  are 
all  in  near  exact  agreement  with  experimental  data^s*®®.  Using  this  surface,  we 
have  investigated  the  gas-phase  bimolecular  reaction  rates  and  mechanism  for 
the  reaction  of  F2  with  C2H4,  reactions  R26-R28®®,  as  well  as  rates  for  the 
unimolecular  decomposition  of  1 ,2-difluoroethane  upon  random  internal 
excitationfo®,  and  upon  mode-specific  excitation^Of. 

There  are  numerous  benchmarks  that  illustrate  the  accuracy  of  our 
results.  Our  calculations  on  gas-phase  C2H4  +  F2  reactions®®  predict  that  the 
reaction  products  are,  in  order  of  importance,  fluoroethyl  radicals  +  fluorine 
atoms  and  HF  +  fluoroethylene.  This  is  in  complete  accord  with  the 


experimental  results  reported  by  Kapralova  et  al.^^  and  with  the  crossed 
molecular  beam  data  obtained  by  Grover  ef  on  the  analogous  benzene  + 
F2  reaction.  For  unimolecular  dissociation  of  1 ,2-difluoroethane,  the 
decomposition  channels  are,  in  order  of  importance,  reactions  R29,  R32,  and 
R31.  The  calculated  microcanonical  rate  coefficient  for  HF  elimination  (R29)  at 
138  kcai/mol  internal  excitation  is  2.5  x  10^^  s*T  The  extrapolated  experimental 
result  reported  by  Chang  et  is  4.0  x  10''^  s-1.  The  calculations  predict 
that  71-74%  of  the  reaction  exothermicity  will  remain  in  the  olefin  whenever 
reaction  R29  occurs.  This  is  in  general  agreement  with  severai  experimental 
studiesi03-i05  Table  I  shows  a  comparison  of  the  relative  population  densities 
of  the  different  HF  vibrational  states  predicted  by  the  calculationsioo-^O''  with 
very  recent  data  reported  by  Arunan,  Wategaonkar,  and  Setserf06  for  the 
analogous  reaction  of  CH3-CF3.  The  extent  of  agreement  is  clearly  very  good 
considering  the  fact  that  the  molecules  involved,  although  very  similar,  are  not 
identical. 

Table  I 

Comparison  of  relative  population  density  of  different 
vibrational  states  of  HF  after  four-center  elimination  reaction 


HF(v)  Cal. 

[CH2F-CH2F]  (a) 

Expt.  [CH3-CF3]  (b) 

0 

2.05 

1.42 

1 

1.00 

1.00 

2 

0.71 

0.60 

3 

0.19 

0.23 

4 

0.05 

0.06 

5 

0.00 

0.00 

(a)  Reference  1 00  (b)  Reference  1 06 


Summary: 

A  global  potential-energy  surface  has  been  developed  for  the 
C2H4F2  system^o,2i  g^d  reactions  of  F2  with  ethylene^^  and  of  the 
unimolecular  decomposition  of  1,2-dlfluoroethane^^^’^o^  have 
been  investigated  to  provide  a  means  of  assessing  the  accuracy  of 
our  computation  methods  for  systems  of  this  degree  of  complexity 
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and  to  provide  the  necessary  background  for  determination  of  the 
effects  of  dose  confinement  upon  chemicai  reactions,  which  is  to 
be  investigated  in  future  studies.  Comparison  of  the  results  with  a 
wide  variety  of  data  obtained  from  thermal,  molecular  beam,  and 
chemical  activation  experiments  shows  the  accuracy  of  the 
calculations  to  be  sufficient  to  draw  meaningful  conclusions  about 
the  dynamics. 

iL£. _ Intramolecular  Energy  Transfer  Rates: 

In  many  theoretical  investigations,  it  is  necessary  to  compute  the  rates 
and  pathways  of  intramolecular  energy  flow.  Such  calculations  have  proven  to 
be  difficult,  especially  in  highly  coupled,  polyatomic  systems.  In  general, 
classical  studies  of  intramolecular  energy  transfer  involve  the  integration  of  the 
Hamiltonian  equations  of  motion  on  some  potential-energy  surface.  "Bond"  or 
"mode  energies"  are  then  computed  from  the  results.  Intramolecular  energy 
transfer  pathways  and  rates  are  inferred  from  the  calculated  time  variation  of 
these  quantities ^0^.  This  procedure  obviously  involves  an  arbitrary  definition  of 
the  "bond  energy"  which  generally  assumes  a  mode  separability  that  does  not 
exist.  Consequently,  all  potential  and  kinetic  coupling  terms  involving  the  mode 
coordinates  are  omitted  from  the  definition.  As  a  result,  one  can  never  be 
certain  whether  a  variation  in  "bond  energy"  is  due  to  actual  energy  transfer  to 
or  from  other  modes  or  merely  to  changes  in  the  magnitudes  of  the  omitted 
coupling  terms.  Nor  is  it  possible  to  be  certain  that  the  results  themselves  are 
not  dependent  upon  the  arbitrary  definition  adopted  for  the  "bond"  or  "mode 
energy". 

We  have  recently  reported  a  general  method  for  analyzing  the  results  of 
classical  trajectory  calculations  to  obtain  the  details  of  intramolecular  energy 
transfer  that  obviates  the  need  to  arbitrarily  define  a  "bond"  or  "mode  energy"io. 
The  method  is  based  on  the  determination  of  the  time  dependence  of  the 
normal-mode  velocities  by  projection  of  the  instantaneous  Cartesian  velocities 
of  the  atoms  onto  the  normal-mode  vectors.  The  basic  idea  is  as  follows: 

Let  Lj  (i  =  1 ,2,3,...,3N)  represent  a  set  of  normalized  (3N  x  1) 
transformation  vectors  that  project  the  normal-mode  vibrations  (1  <  i  <  3N-6), 
the  center-of-mass  translations  (3N-5  <  i  <  3N-3),  and  the  rotations  about  the 
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molecular  center-of-mass  (3N-2  <  i  <  3N)  onto  the  Cartesian  displacements,  q\ 
(j  =  1 ,2,3,...,3N).  That  is,  at  time  t,  the  instantaneous  Cartesian  velocities  may  be 
written  as  linear  combinations  of  the  corresponding  elements  of  the  Lj : 

(18)  qj(t)  =  ZQj(t)Lij  for(i  =  1,2,3.....3N). 

i 

Equation  (18)  may  be  written  in  matrix  form  as 

(19)  q(t)  =  LQ(t), 

where  L  is  a  (3N  x  3N)  square  matrix  whose  columns  are  the  normalized 
projection  vectors  Lj.  q(t)  is  a  (3N  x  1)  column  vector  whose  elements  are  the 
Cartesian  velocities  and  Q(t)  is  a  (3N  x  1)  column  vector  whose  elements  are 
the  normal  mode,  center-of-mass,  and  rotation  velocities.  The  normal-mode 
velocities  may  therefore  be  obtained  from 

(20)  Q  (t)  =  L-'  q  (t)  . 

The  kinetic  energy,  T,  at  time  t  is 

(21)  T(t)  =  0.5Smi42(t). 

Substitution  of  Eq.(18)  yieids 

(22)  T  (t)  =  0.5  [  2  mi  1 1  Qj  (t)  Qr  (t)  Lij  U  ]  =  0.5  [  Z  mi  I  Qj2  (t)  Lij2 1 

i  j  k  i  j 

since  the  kinetic  energy  is  diagonal  when  expressed  in  terms  of  the  normal 
mode  velocities.  Consequently,  the  kinetic  energy  may  be  written  as 

(23)  T(t)  =ZakQk2(t), 

k 

provided  the  mode  constant  an  is  defined  to  be 

(24)  an  =  0.5  Z  mi  Lik2  . 

i 

Equation  (23)  shows  that  the  total  molecular  kinetic  energy  is  the  uncoupled 
sum  of  the  individual  mode  kinetic  energies,  Ti  (t),  where 

(25)  Ti(t)  =  aiQi2(t)  . 

Since  the  total  energy  associated  with  a  given  mode  alternates  between 
potential  and  kinetic  energy  with  a  frequency  closely  associated  with  the  mode 
frequency,  the  flow  of  energy  through  the  molecule  may  be  conveniently  and 
exactly  monitored  by  examination  of  the  envelope  functions  for  plots  of  the  time 
variation  of  the  various  Tj(t).  Thus,  for  the  special  case  of  an  uncoupled 
harmonic  potential  where  the  normal  modes  energies  are  constants  of  the 
motion,  we  would  expect  to  see  Ti(t)  oscillate  with  a  frequency  characteristic  of 
the  mode  frequency  but  with  an  amplitude  that  is  independent  of  time.  The 


envelope  function  would  be  a  constant  reflecting  the  constant  energy  in  the 
mode.  If,  on  the  other  hand,  energy  is  rapidly  leaving  a  particular  mode,  this  will 
be  reflected  by  an  envelope  function  with  a  large  negative  slope.  This  will  be 
the  case  independent  of  the  nature  of  the  potential-energy  surface. 

The  average  total  energy  associated  with  a  given  mode  may  be  obtained 
from  the  results  of  the  trajectory  integration  using  Eqs.(20)  and  (25)  in 
conjunction  with  the  virial  theorem  which  requires  that  the  total  energy,  on  the 
average,  be  equipartitioned  between  kinetic  and  potential  energy.  For  a  given 
potential,  the  results  are  exact  within  the  framework  of  the  classical 
approximation. 

This  method  has  been  applied  to  the  study  of  intramolecular  energy  flow 
in  gas-phase  1 ,2-difluoroethaneio.  Decay  rates  and  pathways  of  energy  flow 
for  initial  excitation  of  each  of  the  18  vibrational  modes  have  been  determined. 
The  results  obtained  from  the  time  variation  of  the  normal-mode  velocities  were 
used  to  extract  the  entire  first-order  mode-to-mode  energy  transfer  rate 
coefficient  matrix.  Essentially,  this  provides  a  complete  description  of 
intramolecular  energy  flow  in  the  molecule. 

Nummary; 

We  have  developed  an  exact  projection  method  that  permits 
intramolecular  energy  transfer  rates  and  pathways  to  be 
determined  from  the  results  of  classical  trajectory  calculations 
without  the  necessity  to  arbitrarily  define  a  "bond"  or  "mode 
energy".  The  method  Is  based  upon  the  calculation  of  the  envelope 
functions  of  the  mode  kinetic  energies  in  a  coordinate  system  in 
which  the  kinetic  energy  is  diagonal.  The  method  has  been  applied 
to  the  determination  of  IVR  rates  and  pathways  in  1,2- 
difluoroethane  ^  o. 
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